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Nuclear reactors are one of the most intense, pure, controllable, cost-effective, and well-understood 
sources of neutrinos. Reactors have played a major role in the study of neutrino oscillations, a phe¬ 
nomenon that indicates that neutrinos have mass and that neutrino flavors are quantum mechanical 
mixtures. Over the past several decades reactors were used in the discovery of neutrinos, were 
crucial in solving the solar neutrino puzzle, and allowed the determination of the smallest mixing 
angle O 13 . In the near future, reactors will help to determine the neutrino mass hierarchy and to 
solve the puzzling issue of sterile neutrinos. 


I. INTRODUCTION 

Neutrinos, the products of radioactive decay among 
other things, are somewhat enigmatic, since they can 
travel enormous distances through matter without in¬ 
teracting even once. Understanding their properties in 
detail is fundamentally important. Notwithstanding that 
they are so very difficult to observe, great progress in this 
field has been achieved in recent decades. The study of 
neutrinos is opening a path for the generalization of the 
so-called Standard Model that explains most of what we 
know about elementary particles and their interactions, 
but in the view of most physicists is incomplete. 

The Standard Model of electroweak interactions, de¬ 
veloped in late 1960s, incorporates neutrinos (i^e, 
as left-handed partners of the three families of charged 
leptons (e“, /i“, t~). Since weak interactions are the 
only way neutrinos interact with anything, the unneeded 
right-handed components of the neutrino field are absent 
in the Model by definition and neutrinos are assumed 
to be massless, with the individual lepton number (i.e. 
the number of leptons of a given flavor or family) being 
strictly conserved. This assignment was supported by 
the lack of observation of decays like /r'*' —>■ e'*' -|- 7 or 
Kl ^ + fjN, despite the long search for them. 

The discovery of neutrino oscillations over the past 
several decades proved that these assumptions were in¬ 
correct. That discovery represents one of the very few 
instances that show that the Standard Model is indeed 
incomplete. The phenomenon of neutrino oscillations 
means that neutrinos have a finite mass, albeit very 
small, and that lepton flavor is not a conserved quan¬ 
tity. Box 1 explains the basic physics of neutrino oscilla¬ 
tions and their relation with neutrino masses, and intro¬ 
duces the parameters used in the oscillation formalism. 
Determination of all their values, with ever increasing 
accuracy, was and continues to be the main goal of neu¬ 
trino experiments. The current experimental values of 
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FIG. 1. Illustration of neutrino oscillations. The ex¬ 
pected flavor composition of the reactor neutrino flux, for 
neutrinos of 4 MeV energy used as an example, is plotted as 
a function of distance to the reactor cores. The fraction of 
neutrino flavors is calculated based on the neutrino oscilla¬ 
tion theory introduced in Box 1. Reactor neutrino oscilla¬ 
tion experiments are placed at different baselines to measure 
the oscillation features driven by different mechanisms. The 
experiments are usually categorized as: very short-baseline 
{L ~ 10 m); short-baseline {L ~ 100 m); kilometer-baseline 
(L ~ 1 km); medium-baseline {L ~ 50 km); and long-baseline 
{L > 100 km) experiments. Details of the past and future re¬ 
actor experiments, their goals and achievements, are the main 
topics of this Review. 


the mass-squared differences and the mixing an¬ 

gles 9ij can be found in the latest editions of the Review 
of Particle Physics [1] and is also shown in Box 1. Histor¬ 
ically, the concept of neutrino oscillations was first con¬ 
sidered by Pontecorvo [2, 3] and by Maki, Nakagawa and 
Sakata [4], hence the neutrino mixing matrix is usually 
called the PMNS matrix. 

The study of reactor neutrinos played a very signifi¬ 
cant part in the discovery and detailed study of neutrino 
oscillations and will continue to be essential to its further 
progress. Here we briefly review the main points of this 
saga. Fig. 1 illustrates how the flavor composition of the 
reactor neutrino flux, starting as pure at production 
(see the next section for details), is expected to oscillate 
as a function of distance. Experimental verification of 
this behavior, and quantitative analysis of the results, 
are the main topics discussed below. 
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The existence of neutrinos was predicted by Pauli al¬ 
ready in 1930 [5] in his famous letter attempting to ex¬ 
plain the continuous electron energy distribution in nu¬ 
clear beta decay. However, it was not until 1953-1959 
that Reines and Cowan [6-8] were able to show that neu¬ 
trinos were real particles. Their observation used the 
electron antineutrinos emitted by a nuclear reactor and 
started a long tradition of fundamental discoveries using 
reactor Pe’s. 

In the early experiments detectors were placed at dis¬ 
tances L < 100 m [9-15] (for a review see [16]). These 
pioneering short-baseline experiments, in agreement with 
the later established three-neutrino oscillation theory, did 
not observe variations with distance, but they were im¬ 
portant for the understanding of the reactor neutrino flux 
and spectrum. The KamLAND experiment [17-19] in 
the 2000s convincingly showed that the earlier solar neu¬ 
trino measurements were indeed caused by oscillations. It 
demonstrated for the hrst time that the reactor neutri¬ 
nos indeed oscillate, i.e. that the i>e component changes 
with LjE^^ as explained in Box 1. It also allowed the 
most accurate determination of the mass-squared differ¬ 
ence Am^i. 

In the next generation of reactor experiments including 
Daya Bay [20, 21], RENO [22] and Double Chooz [23, 24], 
the longstanding puzzle of the value of the mixing angle 
013 was successfully resolved; it turns out that its value 
013 ~ 8.9° (or sin^ 20i3 = 0.093 [1]) is not as small as 
many physicists expected. That discovery opened op¬ 
portunities for further experiments that should eventu¬ 
ally let us determine the so-far missing fundamental fea¬ 
tures of the oscillations, the neutrino mass hierarchy, and 
the phase 5cp that characterizes the possible CP (charge 
and parity) violation. The planned reactor experiments, 
JUNO [25] and RENO-50 [26], promise to be an impor¬ 
tant step on this path. 

Most of the oscillation results are well described by 
the simple three-neutrino generation hypothesis. How¬ 
ever, there are a few anomalous indications, the so-called 
reactor antineutrino anomaly [27] among them, that can¬ 
not be explained this way. If conhrmed, they would in¬ 
dicate the existence of additional fourth, or fifth, and so 
on neutrino families called sterile neutrinos. These neu¬ 
trinos lack weak interactions and would be observable 
only when mixing with the familiar active neutrinos. The 
proposed very short-baseline reactor experiments at dis¬ 
tances ^10 m will test whether this fascinating possibility 
is realistic or not. 

The discovery of neutrino oscillations is one of the most 
important events in the field of particle physics at the 
present time. In this work we briefly review the contri¬ 
bution to this achievement of the experiments involving 
neutrinos emitted by nuclear reactors. First, we show 
how the reactor flux and its spectrum are determined. 
Then, we describe the success of the KamLAND exper¬ 
iment that is complementary to the exploration of solar 
neutrino oscillations, and the determination of the small¬ 
est mixing angle 0 i 3 using three independent reactor ex¬ 


periments. 

While those achievements involved considerable effort, 
answering the remaining open questions is even more 
complicated. We begin by describing the planned large 
reactor experiments, JUNO and RENO-50, aiming at the 
difficult determination of the so far unknown neutrino 
mass hierarchy or mass ordering. We then touch on the 
future very short-baseline reactor experiments aiming to 
test the tentative and unexpected possibility that addi¬ 
tional light sterile neutrinos might exist. 

II. REACTOR NEUTRINO FLUX AND 
SPECTRUM 

Nuclear reactors derive their power from fission. The 
fission fragments are neutron rich and undergo a cascade 
of (3 decays. Each fission is accompanied by approxi¬ 
mately 6 j5 decays, producing an electron and electron 
antineutrino each. The decay energy, typical for the nu¬ 
clear j3 decay, is a few MeVs, rarely exceeding ^8 MeV. 
Since a typical power reactor core has thermal power 
of about 3 GWt/i, and produces ~ 200 MeV of energy 
in each fission, the typical yield of Ve at equilibrium is 
~ 6 X 10^°i>e core“^ s“^. Reactors are therefore powerful 
sources of low energy Pe’s. 

Neutrinos can be detected through charged current in¬ 
teractions when they produce charged particles, electron 
(mass 0.511 MeV), muon (mass 105.7 MeV) or tan (mass 
1776.8 MeV), with neutrino energy sufficient to produce 
them. The reactor energy is low, thus only reactions 
producing positrons are possible. Hence, to study neu¬ 
trino oscillations with nuclear reactors, one must use the 
disappearance type of tests, i.e. measure the flux as a 
function of the distance L and energy (see Box 1 for 
the detailed formalism) and look for the deviation from 
the simple geometrical scaling. Traditionally such mea¬ 
surements were compared with the expected Ve spectrum 
of the reactor. Good knowledge of that spectrum, its nor¬ 
malization and the associated uncertainties, is essential 
in that case. To reduce the dependence on the knowl¬ 
edge of the reactor spectrum, more recent experiments 
[20, 22] use two essentially similar detectors, one nearer 
the reactor and another farther away. 

There are two principal and complementary ways to 
evaluate antineutrino spectra associated with fission. 
The summation method uses known cumulative fission 
yields V„(Z, A,t), and combines them with the experi¬ 
mentally known (or theoretically deduced) branching ra¬ 
tios bn^i{EQ) of all /3-decay branches with endpoints Eq 
and a normalized shape function of each of these many 
thousands of /3 decays, Pjj{Ejy,EQ, Z), 

dN 

— = E„r„(Z,A,t)E,6„,,(£;*)Pp(Up,U*,Z) . (1) 

There are several difficulties with this method. The 
branching ratios and endpoint energies are sometimes 
poorly known at best, in particular for the short-lived 
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FIG. 2. Detection of reactor v^. In the bottom of the 
figure, the reactor he flux from the individual isotopes [28, 29], 
weighted by their typical contribution to the total flux in a 
commercial reactor, is shown. The detection of he relies on 
the inverse beta decay reaction, whose cross section [30, 31] 
is shown as the blue curve. Their product is the interaction 
spectrum measured by the detectors, shown as the red curve. 
The steps involved in the detection are schematically drawn in 
the top of the figure. The he interacts with a proton, becoming 
a positron (e”*") and a neutron. The e”*" quickly deposits its 
energy and annihilates into two 511-keV y-rays, which gives 
a prompt signal. The neutron scatters in the detector until 
being thermalized. It is then captured by a proton ~ 200 
/is later and releases a 2.2-MeV y-ray (the capture time can 
be significantly reduced by the doping of isotopes with very 
large neutron capture cross section such as gadolinium). The 
detection of this prompt-delayed signal pair indicates an he 
candidate. 

fragments with large Q values and many branches. The 
individual spectrum shape functions Pp{Ejj, Eq, Z) re¬ 
quire description of the Coulomb distortions including 
the nuclear finite size effects, weak magnetism, and ra¬ 
diative corrections. In addition, not all decays are of the 
allowed type. There are numerous (about 25%) first for¬ 
bidden decays involving parity change, where the individ¬ 
ual spectrum shapes are much more difficult to evaluate. 

The other method uses the experimentally determined 
spectrum of electrons associated with fission of the prin¬ 
cipal reactor fuels. That spectrum has been measured 
at ILL Grenoble for the thermal neutron fission of 
^^®Pu and ^"^^Pu [32-34] and recently also for the fast 
neutron fission of Munich [35]. These electron 

spectra are then transformed into the Vg spectra using 
the obvious fact that these two leptons share the total 
energy of each /3-decay branch. The transformation is on 
the basis of fitting first the electron spectra to a set of 
30 or more virtual branches, with the equidistant end¬ 
point spacing, determining from the fit their branching 


ratios. The conversion to the Dg spectrum is performed in 
each of these virtual branches. That conversion is based 
on the assumption that the electron spectrum is known 
precisely. When all virtual branches are put together 
one has to also take into account that different nuclear 
charges Z contribute with different weights to different 
electron and Pg energies. While the conversion would in¬ 
troduce only minimum uncertainty if all decays would be 
of the allowed shape, the presence of the first forbidden 
decays introduces additional uncertainty whose magni¬ 
tude is difficult to determine accurately. 

The summation method was used initially in [36-40] 
and in the more recent version in [29]. The conversion 
method was first used in [32-34], more details can be 
found in [41] and the more recent version in [28]. Nat¬ 
urally, the thermal power of the reactor and its time¬ 
changing fuel composition must be known, as must the 
energy associated with fissions of the isotopes 
239pu, 24ipy 238U addition, as already men¬ 
tioned, small corrections to the spectrum shape of indi¬ 
vidual /3-decay branches due to the radiative correction, 
weak magnetism, nuclear size, and so on must be cor¬ 
rectly included. Difficult to do accurately, but of a par¬ 
ticular importance, is to take into account the spectrum 
shape of the numerous first forbidden /3 decays [42]. The 
overall uncertainty in the flux was estimated in [28, 29] 
to be ^ 2%, but when the first forbidden decays are in¬ 
cluded it is estimated in Ref. [42] that the uncertainty 
increases to ~ 5%. 

In essentially all reactor neutrino oscillation studies, 
the Vg are detected using the inverse neutron /3-decay 
reaction 

Dg +p e'^ +n , a = 9.53^'^^(1 -bcorr.) x 10“^^cm^ , 

( 2 ) 

whose cross section is accurately known [30, 31] and de¬ 
pends primarily on the known neutron decay half-life. 
(At the same time the recoil, radiative corrections etc., 
must be also be taken into account.) Since the neu¬ 
tron is so much heavier than the available energy, its 
kinetic energy is quite small (tens of keV) and thus the 
principal observables are the number and energy of the 
positrons. Most importantly, the correlated observation 
of the positrons and the delayed neutron captures is a 
powerful tool for background suppression. Note that the 
reaction (2) has a threshold of 1.8 MeV, only Pg with 
energy larger than that can produce positrons. 

In Fig. 2 we illustrate the energy dependence of the 
reactor Pg flux, the detection reaction cross section and 
their product, i.e. the measured antineutrino spectrum. 
The contributions of the individual isotopes to the Pg 
flux, weighted by their typical contribution to the reac¬ 
tor power are also shown. The top part of the figure 
schematically indicates the steps involved in the Pg cap¬ 
ture on proton reaction. 
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Box 1: Neutrino Oscillations 

Neutrinos are produced with a definite ‘fiavor’: or v^-. For example, in nuclear /3 decay electron 

antineutrinos (i>e) are always produced together with an electron. Similarly, if a positively charged muon (/r'*') 
is produced in the decay of the meson tt"*', the muon neutrino (y^) is always produced as well. However, if 
neutrinos have a finite mass, the fiavor composition of a neutrino beam could vary regularly as a function of the 
distance and energy. This behavior, called neutrino oscillation, is a subtle consequence of quantum mechanics 
that postulates that a neutrino of a given fiavor need not be a state of a definite mass, but instead could be a 
coherent superposition of several states of definite masses. Here we explain the basic ideas of this phenomenon. 

Let us assume that there are only two massive neutrinos Viy = 1,2 with different masses m^. In quantum 
mechanics the development of a state ip with momentum p is governed by il){t) = (We use the usual 

system of units where h = c = 1). For highly relativistic neutrinos p = y/— m? E — w?/2E. Thus, 
when neutrinos propagate in vacuum over a distance L, each acquires the phase i'i{L) = Vi^O) exp{—imfL/2E). 
(The overall phase is skipped. Additional phases are acquired when neutrinos propagate in matter, the so-called 
“MSW effect” [43, 44], which for simplicity will not be discussed here.) Assume further that the fiavor neutrinos 
Vg and Va, i.e. the neutrinos that are the partners of charged leptons in the weak interactions, are coherent 
superpositions of the states Vi, i.e. Vg = cos 0 i^i -|- sin 0 j^ 2 , and analogous but orthogonal combination represents 
the other fiavor neutrino Va = —sinOvi -\- cos9v2- These mixtures are characterized by the parameter 0, the 
so-called mixing angle. 

Consider now a beam of neutrinos that at L = 0 is pure Then 

Ve{L) = cos0e-*™i^/2^:/i(O) + sin0e-™=^/2^z/2(O) . (3) 

In order to observe this beam at L we need to use the weak interactions. We must therefore project the Vi back 
to the flavor basis Ve and Thus 

Ve{L) = [cos2 6le-™N/2i5 + sin2 6le-™2i'/2i5]^^(0)-sin6lcos6l[e-*’"i^/2-®-e-*'"=^/2i5]j,^(0) . (4) 

The probability that we detect Vf, at the distance L is just the square of the corresponding iZe(O) amplitude. The 
probability of detecting is the square of the J^ct(O) amplitude. After some simple algebra this becomes: 

P{ve —>• i/g) = 1 — sin^ 20sin^ , P{ve —> Va) = sin^ 20sin^ , (5) 

4ii/ 4A 

where — m\ is the difference of the squares of the neutrino masses. 

We see that, provided Am^ ^ 0 and 0 7 ^ 0 or 7 r/ 2 , the composition of the neutrino beam oscillates as a function 
of LjEl, with the amplitude sin^ 20 and wavelength 

r ^ _ r 2.48A(MeV) 

Am 2 (eV 2 ) ' 

Observation of neutrino oscillations, consequently, constitutes a proof that at least some of the neutrinos have 
a finite mass and that the superposition is a nontrivial one. Generalization to the realistic case of three neutrino 
flavors and three states of definite mass is straightforward. The corresponding mixing is then characterized by 
three mixing angles ^ 12 , 6 * 13 , ^* 23 , one possible CP violating phase Scp^ and two mass square differences Am^i 
and Am| 2 - As of 2014, in the latest edition of the Review of Particle Physics [1], the best measured values 
using existing data are sin^( 20 i 2 ) = 0.846 ± 0.021, sin^( 20 i 3 ) = 0.093 ± 0.008, sin^(2023) = O.QQOIq;™^, Ato^i = 
(7.53 ±0.18) X 10“®eV^, and Am^2 = (2.44 ±0.06) x 10“^eV^ (assuming m3 > m2)- The mass ordering between 
m 3 and m 2 , often referred to as the neutrino mass hierarchy, and the value of Scp, are currently still unknown. 

Neutrino oscillations have no classical analogue. They are purely quantum-mechanical phenomena, the conse¬ 
quence of the coherence of the superposition of states. 
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III. EXPLORING SOLAR NEUTRINO 
OSCILLATIONS ON EARTH 

Since the late 1960s, a series of solar neutrino experi¬ 
ments [45-49], using charged current reactions, have ob¬ 
served a large deficit of solar z/g flux relative to the Stan¬ 
dard Solar Model (SSM) [50] prediction. It appeared that 
more than half of the solar neutrinos were missing. This 
was referred to as the “Solar Neutrino Problem”. In 2001, 
the SNO solar neutrino experiment [51], for the first time, 
successfully measured the total flux of all three neutrino 
flavors {ve, and Vr) through the neutral current chan¬ 
nel v + d —>■ V -\-p + n using heavy water as a target, that 
yielded results consistent with the SSM. The SNO result 
is considered the “smoking gun” evidence of the neutrino 
oscillation explanation to the Solar Neutrino Problem — 
the solar neutrinos, produced as electron-neutrinos from 
fusion and other reactions in the central region of the 
Sun, are transformed into other flavors when they arrive 
at the Earth. 

The solar neutrino experiments allowed several pos¬ 
sible solutions in the oscillation parameter space of 0 i 2 
and Ato 2 i- A precise measurement of these parameters 
and demonstration of the oscillatory feature, however, 
was hindered by the relatively large uncertainty of the 
solar z/e flux predicted by the SSM, the large matter ef¬ 
fect inside the Sun, and the extremely long distance the 
neutrinos travel. A reactor neutrino experiment, mea¬ 
suring the same disappearance channel as the solar neu¬ 
trino experiments assuming CPT invariance, overcomes 
these difficulties. With well-understood and controllable 
~MeV z/e source, a reactor experiment at ~100 km base¬ 
lines can explore with high precision the so-called “Large 
Mixing Angle (LMA)” parameter region suggested by the 
solar neutrino experiments. To do that, the KamLAND 
experiment [17] was built in early 2000s to explore the 
solar neutrino oscillations on Earth. 

To shield against cosmic rays, the KamLAND detec¬ 
tor was placed at the site of the former Kamiokande ex¬ 
periment [48] under the summit of Mt. Ikenoyama in the 
Japanese Alps. The vertical overburden is 2700-m water- 
equivalent (m.w.e). It is surrounded by 55 Japanese nu¬ 
clear reactor cores, which then produced about 30% of 
the total electricity in Japan. The De flux-weighted aver¬ 
age baseline is about 180 km, well suited for KamLAND 
to study the parameters suggested by the solar neutrino 
experiments. The reactor operation information such as 
thermal power, fuel burn-up, and fuel exchange and en¬ 
richment records were provided by all Japanese reactors, 
which allowed KamLAND to calculate the instantaneous 
fission rate of each isotope accurately. 

The KamLAND detector consists of 1 kton of highly 
purified liquid scintillator (LS), enclosed in a 13-m- 
diameter transparent balloon suspended by ropes in min¬ 
eral oil (MO). The MO is housed inside a 18-m-diameter 
stainless steel (SS) sphere, where an array of 1325 17- 
inch and 554 20-inch photomultiplier tubes (PMTs) is 
mounted. The MO shields the inner LS region from ex¬ 


ternal radiation from PMTs and SS. Purified water (3.2 
kton) is used to provide further shielding against ambi¬ 
ent radiation and operates as an active cosmic muon veto 
detector. With regular central-axis deployments of ra¬ 
dioactive sources and dedicated off-axis deployments [52] , 
KamLA ND achiev ed an excellent position res olution o f 
12 cm/ A/E(MeV), energy resolution of 6.5%/^/^(MeV), 
and absolute energy-scale uncertainty of 1.4%. 

Even with such powerful reactor De sources and a mas¬ 
sive detector, the long baseline suppresses the expected 
signal at KamLAND to only about one reactor z?e event 
per day. In comparison, the background from internal 
and external radioactivity is one million times higher. 
The experiment is only possible owing to the powerful co¬ 
incidence signature (the positron followed by the delayed 
neutron capture 7 ) of inverse beta decay, as illustrated in 
Fig. 2. A time difference of less than 1 millisecond and 
distance less than 2 meters between the prompt and de¬ 
layed signals is required in the analysis. Only the inner¬ 
most 6 -m-radius scintillator region is used to reduce the 
accidental coincidence from external radioactivity. Infor¬ 
mation about the event energy, position, and time were 
used to further reduce the accidental background to ~5% 
of the candidates. 

The other dominant background (~10%) at Kam¬ 
LAND is caused by the ^^C(a,n)^®0 reaction where the 
a-decay comes from ^^°Po, a decay product of ^^^Rn 
that is naturally present in the air and many materi¬ 
als as traces, but is sufficient to induce a measurable 
contamination of the scintillator during its production. 
The neutron scattering off proton or de-excitation 
produces a prompt signal, followed by a neutron capture 
delayed signal. This then mimics a true bg event. The 
rest of the backgrounds include the following: the an¬ 
tineutrinos produced in the decay chains of ^^^Th and 
in the Earth’s mantle and crust, so-called geoneu¬ 
trinos; the cosmogenic beta-delayed neutron emitters ®Li 
and ®He; the fast neutrons from muons passing through 
the surrounding rock, as well as atmospheric neutrinos. 

Fig. 3 (a) shows the prompt energy spectrum of z?e 
candidate events, observed with 2.9 kton-year exposure, 
overlaid with the expected reactor and background 
spectra. A total of 1609 events were observed, which is 
only about 60% of the expected signal if there are no 
oscillations. The ratio of the background-subtracted Ve 
candidate events to no-oscillation expectation is plotted 
in Fig. 3 (b) as a function of L/E,j. The spectrum indi¬ 
cates almost two cycles of the periodic feature expected 
from neutrino oscillations, strongly disfavoring other ex¬ 
planations of the Vg disappearance. 

The KamLAND results [17-19] are highly consistent 
with the solar neutrino experiments, and have pinned 
down the solar neutrino oscillation solution to the LMA 
region. When combined with the results from SNO, 
they yield the most precise measurements of tan^ 612 = 
0.47t[j;[]^ and = 7 . 59 t[]; 2 i x IQ-^ eV^. This is a 

great example of the complementarity between different 
types of experiments. The SNO and KamLAND’s first 
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FIG. 3. Results from KamLAND. (a) The data points show the measured prompt energy spectrum of v^. candidate events. 
The shaded histograms show the expected backgrounds. The expected reactor spectra without oscillation and with best-fit 
oscillation are shown as the dashed histogram. All histograms incorporate the energy-dependent selection efficiency, which is 
shown on the top. Only ~60% of reactor Ve’s are observed relative to the no-oscillation expectation, (b) The data points show 
the ratio of the background-subtracted t'e spectrum to the expectation for no-oscillation as a function of Lq/Ev. Lq is the 
effective baseline taken as a flux-weighted average {Lq = 180 km). The spectrum indicates almost two cycles of periodic feature 
as expected from neutrino oscillations. The oscillation survival probability using the best estimates of 0i2 and |Am|i| is given 
by the blue curve. The curve deviates from the perfect sinusoidal L/E dependence since KamLAND has multiple baselines. 
Figures are reproduced, with permission, from Ref. [19]. 


results came out within approximately 18 months of each 
other, with the solar experiment being more sensitive to 
the mixing angle O 12 and the reactor experiment to the 
mass-squared difference Am^i. The observation of the 
same effect with two different sources on such different 
scales provides compelling evidence for neutrino oscilla¬ 
tions. 


IV. SEARCHING FOR THE SMALLEST 
OSCILLATION ANGLE 

In contrast to the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix in quark mixing, where all three mixing 
angles are very small [ 1 ], the mixing angles in the neu¬ 
trino mixing matrix appear to be large: 623 , measured 
by the atmospheric [53] and long-baseline accelerator [54] 
neutrino experiments, is consistent with 45° which cor¬ 
responds to maximal mixing; and O 12 , measured by the 
solar neutrino experiments and KamLAND, is about 33°. 
It was therefore natural to expect that the third mixing 
angle, ^ 13 , might be of similar magnitude. 

The cleanest way to measure 613 is through kilometer- 
baseline reactor neutrino oscillation experiments. A non¬ 
zero 013 will cause a deficit of Pg flux at approximately 
1-2 kilometer baseline, as indicated in Fig. 1. The size 
of the deficit is directly proportional to the value of 
sin^ 20 i 3 . Unlike accelerator neutrino experiments, the 
reactor measurements are independent of the CP phase 
and 023 ) and only slightly dependent on the neutrino 
mass hierarchy and matter effect. A high precision mea¬ 


surement can therefore be achieved. 

In the 1990s, two first-generation kilometer-baseline re¬ 
actor experiments, CHOOZ [56] and PALO VERDE [57] 
were constructed to measure 0 i 3 . The CHOOZ detector 
was built at a distance of ~1050 m from the two reactors 
of the CHOOZ power plant of Electricite de France in 
the Ardennes region of France. It took data from April 
1997 to July 1998. The PALO VERDE detector was 
built at distances of 750, 890 and 890 m from the three 
reactors of the Palo Verde Nuclear Generating Station 
in the Arizona desert of the United States. It took data 
between October 1998 and July 2000. Surprisingly, nei¬ 
ther experiment was able to observe the Pg deficit caused 
by ^13 oscillation. As a result, only an upper limit of 
sin^ 20 i 3 < 0.10 at 90% C.L. was obtained [56]. 

The null results from CHOOZ and PALO VERDE, 
combined with the measured values of 623 and 0 i 2 , moti¬ 
vated many phenomenological speculations of neutrino 
mixing patterns such as bimaximal and tribimaximal 
mixing [58, 59]. In most of these theories, 0 i 3 is either 
zero or very small. A direct consequence of a vanishing 
013 is that the CP violation in the leptonic sector, even 
if large, can never be observed in the neutrino oscilla¬ 
tion experiments. The importance of knowing the pre¬ 
cise value of 013 provoked a series of world-wide second- 
generation kilometer-baseline reactor experiments in the 
twenty-first century, including Double Chooz [23] in 
France, RENO [22] in Korea and Daya Bay [20] in China, 
to push the sensitivity to 0 i 3 considerably below 10°. Ta¬ 
ble I summarizes some of the key parameters of the five 
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FIG. 4. Results from Daya Bay. (a) The data points show the ratio of the detected to expected no-oscillation Ve. signals at 
the 6 antineutrino detectors located in three experimental halls (EHs) as a function of the weighted baseline. An ~6% signal 
deficit at the far hall relative to the near halls is observed. The oscillation survival probability at the best-fit value is given 
by the red curve. The no-oscillation hypothesis (0i3 = 0) is excluded at 5.2 standard deviations, as shown in the inset, (b) 
The data points show the ratio of the background-subtracted he spectrum to the expectation for no-oscillation in the three 
experimental halls as a function of LeStjE^. The effective baseline Left is determined for each experimental hall (EH) equating 
the multi-core oscillated flux to an effective oscillated flux from a single baseline. A near-complete cycle of the expected periodic 
oscillation feature is observed. The oscillation survival probability using the best estimates of and |Am|i| is given by the 
red curve. Figures are reproduced, with permission, from Ref. [20] and [21]. 



Power 

Baseline Mass 

Overburden 


(GWth) 

(m) 

(ton) 

(m.w.e) 

CHOOZ [56] 

8.5 

1050 

5 

300 

PALO VERDE [57] 

11.6 

750-890 

12 

32 

Double Chooz [23] 

8.5 

400 

8 

120 



1050 

8 

300 

RENO [22] 

16.8 

290 

16 

120 



1380 

16 

450 

Daya Bay [20] 

17.4 

360 

2x20 

250 



500 

2x20 

265 



1580 

4x20 

860 


TABLE I. Key parameters of the reactor Oiz experi¬ 
ments. The table summarizes the key parameters of the five 
past and present reactor 0 i 3 experiments, including the reac¬ 
tor thermal power (in giga-watts), distance to reactors, target 
mass of the detectors, and overburden of the underground site 
(in meter-water-equivalent). 


aforementioned experiments. 

A common technology used in both the first and second 
generation experiments is the gadolinium-loaded liquid 
scintillator as the detection target. Gd has a high ther¬ 
mal neutron capture cross section. With ~0.1% gadolin¬ 
ium loading, the neutron capture time is reduced to ~28 
microseconds from ~200 microseconds for the un-loaded 
scintillator (as used in KamLAND). Furthermore, Gd de¬ 
excitation after the neutron capture releases an 8 -MeV 
gamma-ray cascade, which gives a delayed signal well 
above natural radioactivity (in contrast, neutron capture 
on a proton releases a single 2.2-MeV 7 ). The accidental 


coincidence background is therefore drastically reduced. 

Addition of near detectors at baselines of a few hun¬ 
dred meters is the most significant improvement of the 
second-generation experiments over the previous ones. 
As discussed above, the uncertainty in predicting the re¬ 
actor antineutrino flux is relatively large (2-5%). This 
flux uncertainty, however, can be largely eliminated by 
the relative measurement between near and far detec¬ 
tors. The Double Chooz experiment expands GHOOZ 
by adding a near detector at a distance of ^400 m. The 
installation of that near detector, however, was delayed 
due to civil construction. Double Chooz started taking 
data in May 2011 with only a far detector, and used the 
Bugey4 measurement [14] to normalize the reactor flux. 
The RENO experiment was built near the six reactors of 
the Yonggwang nuclear power plant in Korea. The two 
identical detectors were located at 290 and 1380 m, re¬ 
spectively, from the center of the reactor array. RENO 
started taking data in August 2011. The Daya Bay ex¬ 
periment was built near the six reactors of the Daya Bay 
nuclear power plant in southern China. Daya Bay had 
eight identical antineutrino detectors (ADs). Two ADs 
were placed at ^360 m from the two Daya Bay reactor 
cores. Two ADs were placed at ~500 m from the four 
Ling Ao reactor cores. Moreover, four ADs were placed 
at a far site ^1580 m away from the 6 -reactor complex. 
This modular detector design allows Daya Bay to largely 
remove the correlated detector systematics. Daya Bay 
started taking data in December 2011. 

Compared with the first-generation experiments, the 
second-generation experiments have much larger signal 
statistics by utilizing higher power reactors and larger 
detectors. Among them, Daya Bay has the largest re- 
























actor power (17.4 GWth) and target mass (80 tons at 
the far site,) as shown in Table 1. The underground 
sites are much deeper to allow better shielding from cos- 
mogenic background, in particular compared with the 
case of PALO VERDE. Better chemical recipes of the 
gadolinium-loaded liquid scintillator also improve the 
overall detector performance and long term stability. 

The second-generation reactor experiments were a 
huge success. In 2012, all three experiments. Double 
Chooz, Daya Bay and RENO, reported clear evidence 
of Vg disappearance at ^kilometer baselines after only a 
few month’s running [20, 22, 23]. In particular, Daya Bay 
excluded 6*13 = 0 by 5.2 s.d. with 55 days of data [20]. 
Fig. 4 (a) shows the result of this discovery. The ratio 
of the detected to expected no-oscillation Dg signals at 
the 6 detectors located in the three experimental halls 
is plotted as a function of weighted baseline. The signal 
rate at the far site shows an obvious ^ 6 % deficit with 
respect to the near sites, and fits nicely to the theoretical 
oscillation curve (in red). The precision of the ^13 mea¬ 
surement improved quickly with more data. With the 
data collected in Daya Bay through November 2013 [55], 
the best-fit value is sin^ 2^13 = 0.084 ± 0.005. Although 
the last known, the precision in ^13 measurement ( 6 %) is 
now the best among all three mixing angles. 

Similar to KamLAND, the ratio of the detected Pg 
events to no-oscillation expectation at Daya Bay is plot¬ 
ted in Fig. 4 (b) as a function of LjE^. The combined 
data from the three experimental halls show a near- 
complete cycle of the expected periodic oscillation fea¬ 
ture. The smaller amplitude and shorter wavelength of 
the oscillation, compared to the case of KamLAND, in¬ 
dicate the different oscillation component driven by ^13 
and Am|i. The best-fit frequency of the oscillation yields 
IATO 32 I = 2.391q;]]() X 10“^ eV^ (assuming normal mass 
hierarchy), which is consistent and of comparable pre¬ 
cision with the measurements of accelerator and 
disappearance [60, 61]. By the end of 2017, Daya Bay 
expects to measure both sin^ 2^13 and |Am| 2 l to preci¬ 
sions below 3% [55]. 

The discovery of ^13 represents another good exam¬ 
ple of the complementarity between different types of 
experiments. The first results from the reactor experi¬ 
ments [20, 22, 23] and the accelerator experiments [62, 63] 
were released within approximately 9 months of each 
other, with the reactor experiments measuring Vg dis¬ 
appearance and the accelerator experiments measuring 
Vg appearance. Seeing the same 0 i 3 -driven effects with 
different sources of neutrinos at very different energy and 
baselines is a strong proof of neutrino oscillations. 

The longstanding puzzle of the value of ^13 is now suc¬ 
cessfully resolved. The relatively large value of @13 opens 
the gateway for future experiments to determine the neu¬ 
trino mass hierarchy and to measure the CP-violating 
phase in the leptonic sector. 


V. DETERMINATION OF NEUTRINO MASS 
HIERARCHY 

At present only the absolute values of the neutrino 
mass-squared differences Am |2 and are known, 

not their sign. Depending on whether both Arri^i and 
Ato §2 are positive or whether they are both negative the 
neutrino mass ordering is usually referred to as normal 
or inverted mass hierarchy, respectively. The neutrino 
mass hierarchy (MH) is a problem of fundamental im¬ 
portance [64] that represents an important step in the 
formulation of the Generalized Standard Model of par¬ 
ticle physics. Its determination will reduce the uncer¬ 
tainty in experiments aiming at the measurement of the 
CP-violating phase and it will help in defining the goals 
of the forthcoming neutrinoless double beta decay exper¬ 
iments. It will also improve our understanding of core¬ 
collapse supernovae. 

The reactor Pg-oscillations are modulated by terms 
which depend on Am^i and Am§ 2 - At a medium base¬ 
line of ^60 km, multiple small-amplitude, proportional 
to the sin^ 20 i 3 , oscillation peaks show up on top of the 
long wavelength oscillation with the much larger ampli¬ 
tude proportional to the sin^ 20 i 2 , as shown in Fig. 1. De¬ 
pending on whether the MH is normal or inverted, the 
small-amplitude oscillation pattern shifts slightly. The 
MH information can be extracted from this pattern by 
using a likelihood analysis [65] or the Fourier transform 
method [ 66 , 67]. Additional information regarding the 
neutrino MH could be obtained by combining the re¬ 
actor oscillation analysis with the long-baseline muon 
neutrino disappearance one [ 68 ], as the effective mass- 
squared differences measured there are different combi¬ 
nations of Arri^i, Ato §2 and other oscillation parameters. 

Two medium-baseline reactor experiments, JUNO [25] 
in China and RENO-50 [26] in Korea, have been pro¬ 
posed aiming at the MH determination, among other 
goals. JUNO is currently under construction. The exper¬ 
iment is located in Kaiping city, Guangdong province, in 
southern China. The JUNO detector will be placed un¬ 
derground with a total vertical overburden of 1800 m.w.e. 
JUNO will observe antineutrinos from the Yangjiang nu¬ 
clear power plant (NPP) and the Taishan NPP at equal 
baselines of ^53 km, as illustrated in Fig. 5, near the 
maximal 0 i 2 -oscillation baseline. The Yangjiang NPP 
has six reactors cores of 2.9 GWth each and the Taishan 
NPP has planned four cores of 4.6 CWt^ each, both are 
under construction. The difference between the base¬ 
lines to the two NPPs is controlled to less than 500 m 
to prevent significant degradation of the MH discrimina¬ 
tion power [65]. JUNO is expected to start data taking 
in 2020. The proposed RENO-50 [26] experiment will be 
located in the city of Naju, about 47 km from the Hanbit 
nuclear power plant with six cores of 2.8 GWth each. The 
detector will be placed at underground of Mt. Guemseong 
with an overburden of 900 m.w.e. RENO-50 is expected 
to begin data taking in 2021 . 

The medium-baseline reactor experiments need to have 
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massive detectors, approximately 20 kilotons, in order to 
collect sufficient j/e events in a reasonable timescale (a few 
years). In the following we will primarily use JUNO as an 
example to illustrate the significant challenges in build¬ 
ing such a large experiment. The preliminary design of 
JUNO includes a central detector submerged in a water 
pool with the muon trackers installed on the top of the 
pool. The water pool is equipped with PMTs and acts as 
an active Cherenkov detector for vetoing muons. It also 
provides passive shielding against the natural radioactiv¬ 
ities from the surrounding rock and air. The top track¬ 
ers provide complementary measurements of the cosmic 
muons. The central detector consists of 20 kton of liquid 
scintillator (LS), contained either in a spherical acrylic 
tank supported by stainless steel frames, or a thin trans¬ 
parent balloon contained in a stainless steel vessel. The 
detector looks similar to the one in SNO or KamLAND, 
but is twenty times larger. In order to collect enough 
light, the central detector is viewed by about 18000 20 - 
inch PMTs. The PMTs have implosion containers to mit¬ 
igate the risk of implosion chain reactions. Taking into 
account the mechanical clearance, the PMTs provide a 
near maximal surface coverage of 75%-78%. RENO-50 
has a similar detector design with 18 kt liquid scintillator 
and 15000 20-inch PMTs. 

An energy resolution better than 3%/ yj ^(MeV) is es¬ 
sential for medium-baseline reactor experiments to main¬ 
tain the MH discrimination ability [65]. To achieve that, 
beside keeping a maximal photocathode coverage, addi¬ 
tional technical improvements are necessary. High quan¬ 
tum efficiency (QE ~ 35%) PMTs are necessary in order 
to increase the light detection efficiency. A new type of 
20-inch micro-channel plate (MCP) PMTs is being de¬ 
veloped for JUNO. The light yield and the optical trans¬ 
parency of the LS also need to be improved. Optimizing 
the concentration of scintillation fluors, purification of 
the raw solvent and fluors, and on-line AI 2 O 3 column fil¬ 
tration have been found effective. LS attenuation length 
of more than 30 m is desired. 

Calibration of the absolute energy scale is crucial. 
In particular, three main effects cause non-linear en¬ 
ergy response of a LS detector: scintillator quenching, 
Cherenkov radiation and possible non-linear electronics 
response. If the energy non-linearity correction has large 
uncertainties, particular residual non-linear shapes may 
fake the oscillation pattern with a wrong mass hierar¬ 
chy [69]. The absolute energy scale uncertainty needs to 
be controlled within a few tenths of a percent, which is 
challenging from the experience of KamLAND [52] and 
Daya Bay [55]. The requirement demands a compre¬ 
hensive calibration program for a large detector such as 
JUNO or RENO-50. 

Background control is demanding, in particular due 
to relatively shallow depth of the experimental sites of 
JUNO and RENO-50. The sources of background are 
similar to those of KamLAND. However, the cosmogenic 
®Li and ®He background is significant, due to the much 
higher muon rate. The ®Li and ®He isotopes are mostly 



FIG. 5. JUNO’s layout and expected signal. JUNO is 
located in southern China at an equal baseline of ~53 km from 
the powerful reactors at Yangjiang and Taishan. The shaded 
histograms in the inset show the expected v,, energy spectra 
at JUNO with and without backgrounds after 6 years’ run¬ 
ning (calculated based on the information given in Ref. [25].) 
The multiple-oscillation structure allows the determination of 
neutrino mass hierarchy and precision measurement of the os¬ 
cillation parameters. 


produced by the muons accompanied by large electro¬ 
magnetic or hadronic showers [70]. In KamLAND, if a 
shower muon is tagged, the whole detector is vetoed for 
2 s. Such a veto strategy will lead to a significant signal 
loss at JUNO and RENO-50. Since the lateral distance of 
the isotopes from the parent muon trajectory is approx¬ 
imately exponential [70], a small veto region along the 
muon track can efficiently remove the background with 
minimal loss of signals. Thus, the ability to track the 
shower muons is essential, which demands new develop¬ 
ments in the muon veto system and improvements on the 
simulations and reconstructions. 

JUNO will observe about 60 reactor Vg, events per day. 
The expected energy spectrum is shown in Fig. 5. The 
sensitivity of the mass hierarchy determination at JUNO 
is estimated to exceed 3cr (for the statistical interpre¬ 
tations see [71, 72]) in 6 years [65, 73]. Assuming that 
the effective mass-squared difference measured by the on¬ 
going accelerator experiments can achieve 1.5-1% preci¬ 
sion [74], the MH sensitivity at JUNO can be improved to 
3.7-4. 4(7 [65]. RENO-50 has similar sensitivity reaches. 
In addition to the MH determination, both JUNO and 
RENO-50 have great potentials in the precision measure¬ 
ments of the neutrino oscillation parameters. The exper¬ 
iments expect to measure Am^i, [Aml^j and sin^ 6*12 to 
precisions better than 1%. This offers a major step to¬ 
ward the unitarity test of the neutrino mixing matrix [75] 
and is important to guide the directions of future exper¬ 
iments and theories. 

The next-generation medium-baseline reactor exper¬ 
iments provide a unique opportunity to determine the 
neutrino mass hierarchy with the precision measurement 
of the reactor neutrino spectrum. Most systematic effects 
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are well-understood and studied, although the technical 
challenges are significant. The MH sensitivity is expected 
to reach 3-4 (t. The reactor measurements are indepen¬ 
dent of 023, the CP-violating phase, and the matter ef¬ 
fect. Combining with the future long-baseline accelera¬ 
tor [76, 77] and atmospheric [78, 79] neutrino oscillation 
programs, we will once again have complementary mea¬ 
surements of the neutrino mass hierarchy with different 
types of experiments. Such complementarity has proved 
essential in the history of establishing the phenomenon 
of neutrino oscillations. 


VI. SEARCHING FOR STERILE NEUTRINOS 

Precision electroweak measurements of the decay 
width of the Z boson determine the number of active 
light neutrinos. The result, 2.92 ± 0.05, [80] is obvi¬ 

ously compatible with the three neutrino flavors. The 
three-neutrino framework has been extremely success¬ 
ful in explaining neutrino oscillation results, since only 
two oscillation frequencies, corresponding to the two 
mass-squared differences (Am^i ~ 7.6 x 10“® eV^ and 
~ 2.4 X 10“^ eV^), were observed by the solar, 
atmospheric, accelerator and reactor neutrino oscillation 
experiments. However, in the 1990s, the LSND exper¬ 
iment [81, 82] reported an anomalous event excess in 
the > Pe appearance channel, which could be inter¬ 
preted as an oscillation with the Am^ ~ 1 eV^. Such 
scale is clearly incompatible with the above Ato^i 
A mgj^. Since the LSND result contradicted the three- 
neutrino framework, it is often referred to as the “LSND 
anomaly”. 

The LSND anomaly indicates the existence of addi¬ 
tional fourth or more neutrino families with masses m ^ 
1 eV. Since these additional neutrinos cannot couple to Z 
bosons, they must lack weak interactions and are there¬ 
fore sterile. Sterile neutrinos are observable only through 
their sub-dominant mixing with the familiar active neu¬ 
trinos. The light sterile neutrinos, coincidentally, are 
also among the leading candidates to resolve outstand¬ 
ing puzzles in astrophysics and cosmology [83-86]. On 
the other hand, the light sterile neutrinos are generally 
not “natural” in the theories that extend the neutrino 
Standard Model. For example, the popular type-I see¬ 
saw model [87-90] , which provides an elegant explanation 
of the small neutrino masses and the matter-antimatter 
asymmetry of the universe [91], predicts only heavy ster¬ 
ile neutrinos {m > I0^° eV). If the light sterile neutrinos 
indeed exist, as LSND indicates, they would suggest new 
frontiers in both experimental and theoretical physics. 

The LSND anomaly so-far remains experimentally un¬ 
confirmed, despite many efforts. However, there are sev¬ 
eral hints supporting LSND’s findings, even though none 
are really conclusive. The MiniBooNE experiment, de¬ 
signed at a similar LjE baseline as LSND using acceler¬ 
ator neutrinos, observed event excess in the —>■ and 
P^ —>■ Pf. appearance channels that have been interpreted 


as consistent with LSND [92, 93]. The GALLEX [94] and 
the SAGE [95] solar neutrino experiments, during their 
calibrations using intense neutrino sources (^^Cr, ^"^Ar), 
observed a ~24% event deficit in the disappearance 
channel. This deficit is often referred to as the “Gallium 
anomaly”. Recently, re-evaluations of the reactor P^ flux 
resulted in an increase in the predicted P^ rate [28, 29] 
(see earlier sections for details). Combining the new pre¬ 
dictions with the reactor experimental data at baselines 
between 10-100 m [9-15] suggests a ^4-6% deficit be¬ 
tween the measured and predicted reactor P^. flux, so- 
called “reactor antineutrino anomaly” [27, 96]. These 
experimental anomalies can be interpreted by light sterile 
neutrinos [97], but might also be caused by an imperfect 
knowledge of the theoretical predictions or experimen¬ 
tal systematics. The preferred region (Am^ ^ 1 eV^ and 
sin^ 20 ~ 0.1), however, is in some tension with the limits 
derived from other appearance [98-101] or disappearance 
searches [15, 102-109]. 

There is, therefore, a strong motivation, despite the 
rather confusing present experimental status, to search 
for the light sterile neutrinos. This led to a high pri¬ 
ority world-wide program [110] with many proposed 
next-generation neutrino oscillation experiments. Dif¬ 
ferent technologies will be used, including short-baseline 
accelerator experiments [111-115] with various neu¬ 
trino beams, ®^Cr (^^'^Ce-^'^'^Pr) v^. {P^) source experi¬ 
ments [116-119] near or inside large LS detectors, as well 
as very short-baseline 10 m) reactor (VSBR) P^ ex¬ 
periments. In order to unambiguously resolve the LSND 
anomaly, the oscillation pattern in the LjE space need 
to be observed, as in KamLAND (Fig. 3) and Daya Bay 
(Fig. 4). VSBR experiments provide unique opportu¬ 
nities to do so given the many advantages provided by 
reactors. 

Multiple VSBR experiments have been proposed glob¬ 
ally; in the U.S. (PROSPEGT [120], NuLat [121]), 
France (NUGIFER [122, 123], STEREO [110]), Rus¬ 
sia (DANSS [124], NEUTRINO-4 [125, 126], POSEI¬ 
DON [127]), U.K. (SOLID [128]), and Korea (HA¬ 
NARO [129]). Table H summarizes some of the key pa¬ 
rameters of the proposed VSBR experiments. The oscil¬ 
lation length of the ~1 eV mass-scale sterile neutrinos is 
about 10 meters for reactor PeS, thus all proposed exper¬ 
iments deploy their detectors at distances of 4-20 m from 
the reactor cores. The reactor cores should preferably 
be compact in size to minimize the oscillations inside the 
cores, so most experiments utilize compact research re¬ 
actors with thermal power of tens of mega-watts. Those 
research reactors are typically highly enriched in in 
contrast to the commercial reactors in the nuclear power 
plants. 

Background control is a challenging task in the VSBR 
experiments. The detectors are typically at shallow 
depth (~10 m.w.e.) constrained by the locations of the re¬ 
actor cores. The cosmic-ray related background is there¬ 
fore high. One advantage of using research reactors is 
that they can be turned on or off on demand, which helps 
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Power 

(MWth) 

Baseline 

(m) 

Mass 

(ton) 

Dopant 

Seg. 

PROSPECT [120] 

85 

6-20 

I&IO 

®Li 

Y 

NuLat [121] 

1500 

3-8 

1.0 

i°B, ®Li 

Y 

NUCIFER [122] 

70 

~7 

0.7 

Gd 

N 

STEREO [110] 

57 

~10 

1.8 

Gd 

N 

DANSS [124] 

3000 

9-12 

0.9 

Gd 

Y 

NEUTRINO-4 [125] 

100 

6-12 

1.5 

Gd 

N 

POSEIDON [127] 

100 

5-8 

1.3 

Gd 

N 

SOLID [128] 

45-80 

6.8 

2.9 

Gd, ®Li 

Y 

HANARO [129] 

30 

6 

~1 

Gd 

Y 


TABLE II. Key parameters of the very short-baseline 
reactor experiments. The table summarizes the key pa¬ 
rameters of the proposed very short-baseline reactor experi¬ 
ments, including reactor thermal power (in mega-watts), dis¬ 
tance to reactors, target mass of the detectors, dopant mate¬ 
rial for neutron capture, and whether or not highly segmented 
detectors are planned. 

to measure the non-reactor background. The reactor- 
related backgrounds, such as fast neutrons and high en¬ 
ergy gamma rays, are however more difficult to determine 
as they appear together with the v^, signals. Sufficient ac¬ 
tive veto and passive shielding are necessary. However, 
given the tight space near the reactor cores, they have to 
be carefully designed. 

As shown in Table II, detectors are typically Gd-loaded 
or ®Li-loaded liquid (or solid) scintillators. The Gd-LS 
technology is mature and a good pulse shape discrimi¬ 
nation (PSD) against the neutron background has been 
demonstrated. An advantage of the ®Li-loaded scintilla¬ 
tor is that the delayed neutron capture process ®Li(n, a)t 
produces an a particle and a triton, instead of a 7 -ray. 
This provides a good localization of the delayed signal 
and an additional PSD against the 7 background. Some 
detectors are highly segmented into small cells in order to 
achieve good position resolution and further background 
rejection by using the multi-cell event topologies. There 
are, however, more inactive layers in the segmented de¬ 
tectors so the edge effects have to be accurately simulated 
and measured. It is also more challenging to perform cal¬ 
ibrations and control the relative variations among cells 
for the segmented detectors. For all detectors, sufficient 
light yield is required to precisely measure the reactor 
i>e spectrum and the possible distortions from neutrino 
oscillations. 

Despite the challenges, very short-baseline reactor ex¬ 
periments provide a great opportunity to observe the 


distinctive feature of the light sterile neutrino oscilla¬ 
tions, due to their extended range of energy (1-8 MeV) 
and baselines (5-20 m). The world-wide next-generation 
VSBR experiments, as shown in Table II, are being ac¬ 
tively considered and pursued. Many of them will begin 
taking data [128] in 2015-16. Within a few years’ run¬ 
ning, they expect to cover the parameter region suggested 
by the experimental anomalies with a sensitivity better 
than 5tT and may tell us whether the fascinating possi¬ 
bility of light sterile neutrinos is true or not. 

VII. OUTLOOK 

Over the past ~60 years, nuclear reactors have proven 
to be one of the most powerful tools to study neutrino 
oscillations, the quantum-mechanical phenomenon that 
requires extensions to the Standard Model. Experiments 
at a few kilometers and at a few hundred kilometers from 
the reactor cores have produced some of the most con¬ 
vincing proofs of neutrino oscillations, by observing the 
oscillatory behavior of the reactor UeS in the LjE domain 
during their propagation. Reactor experiments measured 
several key parameters governing the neutrino mixing, in¬ 
cluding 012 , 013, Ato|i and |Am|i|. They are essential 
in establishing the framework of neutrino oscillations. 

Nuclear reactors will continue to help us uncover more 
facts about neutrinos. In the next ~20 years, the up¬ 
coming next-generation reactor experiments will tell us 
what is the neutrino mass hierarchy and whether or not 
light sterile neutrinos exist. The results will have signif¬ 
icant impact on other future programs such as neutri¬ 
noless double-beta decay experiments, long-baseline ac¬ 
celerator experiments, astrophysics and cosmology. Ulti¬ 
mately, they may hold the key to our deeper understand¬ 
ing of fundamental physics and the universe. 
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